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ABSTRACT 

Context. Modeling the variability of the solar spectral irradiance is a key factor for understanding the solar influence on the 
climate of the Earth. 

Aims. As a first step to calculating the solar spectral irradiance variations we reproduce the solar spectrum for the quiet Sun 
over a broad wavelength range with an emphasis on the UV. 

Methods. We introduce the radiative transfer code COSl which calculates solar synthetic spectra under conditions of non-local 
thermodynamic equilibrium (NLTE) . A self-consistent simultaneous solution of the radiative transfer and the statistical equation 
for the level populations guarantees that the correct physics is considered for wavelength regions where the assumption of local 
thermodynamic equilibrium (LTE) breaks down. The new concept of iterated opacity distribution functions (NLTE-ODFs), 
through which all line opacities are included in the NLTE radiative transfer calculation, is presented. 

Results. We show that it is essential to include the line opacities in the radiative transfer to reproduce the solar spectrum in the 
UV. 

Conclusions. Through the implemented scheme of NLTE-ODFs the COSl code is successful in reproducing the spectral energy 
distribution of the quiet Sun. 

Key words. Radiative transfer - Sun: atmosphere - Sun: chromosphere - Sun: UV radiation - Line: profiles - Atomic data 



1. Introduction 

While the total solar irradiance can not be the 
cause of the climate change ove r the past 20 years 



(jLockwood. M. and Frohhch . 2007f ). it is known to influ- 
ence the Earth's pre-industrial climate (Xabitzke . ,2005; 
iKilifarska fc Haiehl . 120051 : IShindell et aII ■ l2003^ . However, 
it has not yet been fully identified how exactly the irra- 
diance variations affect the climate. It has been argued 
that direct forcing through the variation of the total so- 
lar irradiance is too small to produce the observed cli- 
mate changes as, e.g ., observed during the little ice age 



To understand the solar spectrum variations it is es- 
sential to consider all the relevant physical processes in 
the formation of the solar spectrum. This research has al- 
ready been carried out by several authors who use differ- 
ent approaches and focus on various wavelengths ranges. 
Calculations of solar spectr a in ID over broad wavelength 
ranges were carried out by lKuruc3 ([l991, 2005) based on 
LTE radiative transfer calculations with the ATLAS9 and 
ATLAS12 code i n plan e-parallel sy mmetry. Fu r thermo re, 



Vernazza et al 



AU ende Prieto et al.l (l2003allbh 



Fontenla et aP (ligggL EoO?! ). and lAvrett &: Loesed (120081 ) 



(jShindell et al.l . 120011) . One of the prime candidates for 



an indirect effect is the variation of the UV part of the 
solar spectrum, which varies much more than the total 
irradiance. In particular it has already been shown that 
Hi 121.6 nm (Lyman-a) and the Herzberg band around 
200 nm affect ozone and the temperature in the strato- 
sphere (|Rozanov et al.l . l2006l 120021 : lEgorova et all . 12004 ) . 



calculated solar spectra in NLTE and plane-parallel sym- 
metry for the visible and IR. Improved understand- 
ing of the solar spectrum formation has been achieved 
by the d etailed stud y of individual sp ectral lines by 
Avres et al. (2006) and lCollet et all (|2005t) who calculated 



models out of local thermodynamic equilibrium (NLTE) 
Furthermore ' ^^^^ tTTTT— T"— 71 ^o^-.,! 



20021) 



Hubenv & Lites' (1995"), Uitenbroek' (200A 
Fox et al.. (2004) , and Fontcnla et a l. (1999, 20071) 
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improved the synthetic calcul ations by including partia l 
redistribution (PRD). Finally IShort fc Hauschildti (|2005l ) 
computed model spectra of cool stars in NLTE and spher- 



2 



Haberreiter et al.; NLTE calculations of the solar atmosphere 



ical symmetry. Furthermore, 3D NLTE radiative trans- 
fer calculations based on 3D hyd rodynamic sim u lation s 
have been carr i ed ou t, e.g. by lAsplund et al. ( 2000l ). 
Koesterke et al. (l2008h . Here we limit ourselves to ra- 
diative transfer calculations in ID, as our reconstruc- 
tion approach of the spectral solar irradiance is based on 
ID model atmosphere structures of solar surface features 
and their di s tribution on the solar surface dWenzler et al. , 
2006[ I2OO5I : iHaberreiter et al.L l2005t iKrivova fc Solanki . 
200I). 

The formation of the solar spectrum is, in some spec- 
tral ranges, dominated by an immense number of spectral 
lines, e.g. the iron like elements in the UV, also known as 
UV line haze. Due to the characteristics of the solar tem- 
perature structure there are two effects of the line opacity, 
which depend on the height of line formation, that have 
to be distinguished. One effect is the line blocking in the 
photosphere leading to a decrease of intensity. The other 
effect is the excess of intensity due to emission lines in the 
chromosphere. Because of NLTE effects, i.e. the illumi- 
nation from above, the chromospheric emission lines can 
lead to overionization at lower layers. In this paper we do 
not focus on the chromospheric emission lines apart from 
Lyman-a, as they can only be correctly calculated if all 
atomic processes are treated in full NLTE, but discuss the 
line blocking in t he photosphe re. 

As shown by ICollet et al. the lack of photo- 

spheric line blocking leads to an incorrect excess of ion- 
ization. The authors find differences between calculations 
including and excluding line opacity in NLTE calculations 
with a modified version of the statistical equilibrium code 
MULTI (jCarlssonl Il986l ) for cool stars with solar metallic- 
ity. They show that the inclusion of sampled background 
line opacities decreases the mean intensity field at a given 
depth and thereby reduces the radiative ionization rates. 
The authors find that the decrease of ionization in turn 
changes the line strength and leads to a Fe abundance 
that is 0.1-0.15 dex higher than in calculations without 
line blocking. From this it follows that all opacities con- 
tributing to the line blocking have to be accounted for in 
NLTE radiative transfer calculations. 

There have already been a number of approaches 
for including line opacities in the bac kgrou nd contin- 
uum opacity. Anderson fc Athav ( 1989( ) and Anderson 
( 1989f ) first introduced the concept of super levels that 
allows for numerous levels of certain atoms with simi- 
lar energies, which are assumed to be in LTE with re- 
sp ect to each other. This appro ach was later a dopted 
bv lDreizler fc Werneil (fl993 ) and .Hubenv fc Land (fl995.) , 
who used ODFs for the transitions between the su- 
perlevels. Describing line opacities wi th a Mont e- 



Carlo evaluation has beeii appl ie d bv pchmutz et al . 
(ll99lLISchaerer fc Schrmtp (|l994l ). ISchmut j (|l997l ). and 
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Fig. 1. Photoionization cross section for the third level 
of neutral Mg {Mgl 3 in Table fOll implemented in COSI. 
The solid line is the original data from the Opacity Project 
and the crosses are the data interpolated to the frequency 
grid apphed in COSI. 



de Marco fc Schmut d ( 1999( ) in the calculation of the mass 
loss of hot stars. Furthermore. IShort fc Hauschildt t200jj) 
include more than 100,000 spectral lines or iron-like ions 
that blanket the UV-ban d self-co nsiste ntly in the NLTE 
calculation. Collet et al. ( 2005f l and Avrett fc Loeser 



(j2008l ) include the line opacities in the background con- 
tinuum opacity via sampled line opacities. Our approach 
improves the self-consistency of the NLTE calculation by 
including the line opacities from all known transitions of 
all ions into the radiative transfer calculations by means 
of iterated opacity distribution functions which include 
some NLTE effects, and therefore, we name these func- 
tions NLTE-ODFs. 

The main purpose of this paper is not so much a de- 
tailed analysis of line profiles, but the calculation of re- 
alistic solar energy distribution over a broad wavelength 
range with an emphasis on the UV and to be able to re- 
produce the solar observations and reference spectra over 
a broad wavelength range. 

In the next section we describe the radiative transfer 
code COSI, the model atoms, the model atmosphere, the 
implemented cross sections and the different sets of so- 
lar abundances used in our calculations. Then, in Sect. [3] 
we validate the performance of COSI in LTE against the 
ATLAS12 calculations. In Sect.|3]we describe the new con- 
cept of NLTE-ODFs and present the NLTE calculations 
carried out with COSI. In Sect.lHlthe results are discussed 
and finally in Sect. [6] the conclusions are presented. 

2. COSI 

2.1. Introduction 

COSI (code for Solar Irradiance) is a combination 
of two codes, a mo d el atm osph ere code, developed b y 
Hamann fc Schmutd (|l987t) and ISchmutz etall (|l989l ). 
hereafter COSIMA, and the spectrum synth esis pro- 
gram SYNSPEC, g oing bac k to iHubenvl (Il98ll) and fur- 
ther d eveloped by Hubenv ( 19881 ) and iHubenv fc Land 
(|l995l ). The new method introduced by this paper is 
the use of NLTE-ODFs for including the effects of 
line blanketing replacing the Monte-Carlo evaluation 
of the line opacity that was used in earlier applica- 
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tions dSchmutz et al. . 1991 : Schaerer &: Schmutz . 1994 : 



Schmutd . Il997l : Ide Marco fc Schmutd . Il999l) . COSIMA 
calculates the NLTE populations for a set of specified 
atomic levels by solving the radiative transport equations 
simultaneously with the equations for statistical equilib- 
riu m. The radiative transfer is so lved in spherical symme- 
try (jMihalad . [l97i IPeraiahl |200 ll ) . which yields a more re- 
liable emerging intensity at the limb than a plane-parallel 
geometry and allows t o calculate line of sights at and be- 
yond the solar hmb ( Haberreiter et al.l l2008[ ). the con- 
tribution of which becomes increasingly important for 
UV/EUV wavelengths which are formed higher in the so- 
lar atmosphere. 

The spectral synthesis code takes the NLTE level pop- 
ulations for the explicit levels from COSIMA and assumes 
LTE populations for the non-explicit atomic levels. The 
latter are calculated in LTE relative to the NLTE ground 
states. Population numbers of lines that connect a lower 
explicit NLTE level and a non-explicit upper level are as- 
sumed to be in LTE with respect to the explicit level. 
Partition functions are used for the calculation of the ion- 
ization equilibria in LTE, following the Saha-Equation. 
The partiti on functions fo r the e lements H to Zn they are 
taken from iTraving et aD ( 1966f ) and for the heavier ele- 
ments the coefficients for polynomial fits are taken from 
I Irwin (1981. ) . This allows us to considers all bound-bound 
transitions from hydrogen to thallium, leading to more 
than 2.8 x 10*^ hnes. The hnclist was provided by lKuruc3 
(|200(ih . 



2.2. Atmosphere structures 

For the validation of COSI in LTE we apply the at- 
mosphere structure for the quiet Sun by iKurucz ( 1991 



hereafter ASUN), which is calculated in LTE and has a 
monotonic outward decreasing temperature profile with 
the outer boundary at the temperature minimum. 

The NLTE spectra are calculated with the semi-em- 
pirical structure of the solar atmospher e for avera ge super- 
granule cell interior (Model C) bv iFontenla et al.. (1999), 
further referred to as F1999. 



2.3. Model atoms 

The atomic model of hydrogen consists of negative hydro- 
gen, 10 levels of neutral hydrogen each representing the 
quantum numbers up to n = 10, and the ion. Helium is 
modeled by 10 levels of neutral helium and one level for 
the ground state of singly ionized helium. Furthermore, 
we have implemented all metals up to zinc in the NLTE 
calculation. Currently we account for a total number 61 
ions and 114 explicit NLTE levels as given in Table fATT] 
The metal levels are generally selected according to their 
importance for the solar UV continuum opacity, e.g. their 
ionization edges are in the wavelength range from 100 to 
400 nm. 



All bound-bound and bound-free transitions connect- 
ing the explicit levels are fully included in the NLTE cal- 
culation. In the present implementation the millions of ad- 
ditional transitions are set to LTE relative to the ground 
level of the corresponding ionization stage, which is al- 
lowed to deviate from LTE and explicitly calculated. Thus 
the LTE transitions are to some extent affected by NLTE 
conditions and in turn determine the NLTE-ODFs. 

For line transitions connecting an explicit lower level 
with a non-explicit upper level the latter is calculated in 
LTE relative to the NLTE population of the lower level. 
Here however, an inaccuracy occurs for the chromospheric 
lines that are calculated in LTE, because this leads to an 
overestimation of the line emissivity. To avoid this excess 
of emission in the spectral synthesis the line emissivity for 
each depth point outward of the temperature minimum is 
calculated with the excitation temperature set to the value 
of the temperature minimum, i.e. as e^^x — Bx^Trr^i^i^L.x, 
with B being the Planck function and L the depth point 
index of the atmosphere structure. We are aware of the 
fact that this approach clearly underestimates the chro- 
mospheric emission lines. However, the large excess of the 
emission lines due to their wrong LTE calculation is not 
the reasonable solution either. This shortcoming will be 
addressed in the future by including all the relevant levels 
in the NLTE calculation. 



2.4. Cross sections 
2.4.1. Radiative processes 

The bound-free and free-free radiative processes i nvolv- 
i ng ne gative hydrogen are implemented according to JohnI 
(jl988l ). The photoionization cross section for th e ground 
state of hydrogen is calculated according to iMihalad 
( 1967 ) and the cross sections of th e other hydrogen lev- 
els are adopted from Seaton ( 1960l ). The bound- free ra- 
diative processes for helium are calculated according to 
Koester et al.l (|l985t ) . Free-free cross sections of positively 
charged ions and electrons are calculated with an hydro- 
genic appr oximation using the cro ss sections by iBergeil 
11956) and lKarzas fc Latteil (|l96ll ). 

Table ETTI gives the sources for the photoionization 
cross sections for all the 114 explicit NLTE-levels ac- 
counted for in COSI. For the ions Ci, Nai, Mgi, 
All, Si I, Si and Cai the photoionizati on cross sec- 



tions are taken from the Opacity Project (I Seaton et al 
1994,) a nd the ones for Fei are p rovided by iBautista 
(|l997l) : IBautista fc PradhanI (|l997l ). We linearly inter- 



polate the complex photoionization cross sections to 
the slightly coars e r freq uency grid employed in COSI 
([Haberreiter et al.L 120031 ). As an example, the cross sec- 
tion of Mgi is given in FiglTj As theoretical ionization 
energies are slightly off from the observed ones, we cor- 
rected the theoretical values to agree with the measured 
energies taken from the National Institute of Standard 
and Technology (NIST) Atomic Spectra Database. Also, 
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Fig. 2. Comparison of ATLAS12 calculation by iKurucz ( 20051 ) (dotted line) and COSI calculation (solid line) for 
Hi 410. Inm (panel a) and Nil 676. Snm (panel b), using the same abundances, atmosphere structure, microturbulence 
(1.5 km/s), Gaussian macroturbulcnce (FWHM=1.5 km/s) and rotation broadening (2 km/s). 



Thomson scattering and Rayleigh scattering are included 
in COSI. 



2.4.2. Collisional processes 

In COSI the collisional processes involving negative hy- 
drogen include the collisional detachment by electrons and 
neutral hy drogen and charge neutra lization with protons 
as given bv lLambert fc Pagel ( 1968I) . Eq. (2b,c,e) therein. 
The collisional bound-bound cross sections of neutral hy- 
drogen, the optically permitted transitions of helium, as 
well as of the neutral and ioni z ed me tals are calculated 
using the formula of IJefferied (| 19681 ). p.ll8, Eq. 6.24. 
Collisional processes of neutral hydrogen other than with 
negative hydrogen are not yet included. 



Table 1. List of applied solar abundances relative to H 
by number for the elements (EL) with atomic number A 
and atomic mass mA 



2.5. Solar abundances 

The study of solar abundances is an ongoing research topic 
and various data sets are used by different authors. For the 
syntheti c spectra presen ted in this paper, the abundances 
used by iKurucd |l99ll ) (K91) were applied in the LTE 
COSI calculations using t he atmosphere structu re ASUN. 
The abundances us ed bv iFontenla et al.l (l2007l) (F2007), 
which are based on lAsplund et al.l ( 20051 ). are employed 
for the calculations using the F1999 atmosphere structure. 
Table[l] lists both abundance sets and their ratio for the 
elements implemented in the atomic model. K91 is used in 
the COSI LTE calculations. F2007 is the abundance used 
in all the NLTE calculations presented in this paper. 

3. Code validation in LTE 

In order to validate the overall performance of COSI 
in LTE we calculated the emergent flux of the quiet 
Sun with COSI, applying the same atmosphere struc- 
ture (ASUN), abundances (K91), and values for mi- 
croturbulence (1.5 km/s), macroturbulence (Gaussian 
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FWIIM=1.5 k m/s) and rotation broadening (2 km/s) as 
Kurucd (|2005f ) for the ATLAS12 calculations. The rota- 
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Fig. 3. Panel (a) shows the comparison of the continuum 
calculated with the ATLAS 12 code (dotted line) and con- 
tinuum plus hydrogen lines calculated with COSI in LTE 
(solid line), based on identical atmosphere structure and 
abundances. The difference shortward of the Balmer jump 
at 364.5 nm is mainly due to the photoionization cross sec- 
tion of the third Mg i level, shown in Fig.[l] apparently not 
accounted for in the ATLAS 12 calculation. Panel (b) gives 
the corresponding synthetic LTE spectra averaged with a 
1-nm boxcar. Panel (c) shows the ratio between the COSI 
and the ATLAS 12 calculation. 



tion b roadening has been implemented as given by iGravl 
(|1992'). Fig.dlshows the calculations with ATLAS12 (dot- 
ted hne) and COSI (solid hue) for the spectrum centered 



at Hi 410.1nm (panel a) and NiiA676.8nm (panel b). 
The line profiles of Hi410.1nm and in particular the 
wings of this line are almost identical, showing that in 
COSI the radiative transfer and the Stark broadening are 
computed consistently with the ATLAS12 code. The ab- 
solute flux level of the continuum at the Nil wavelength 
is slightly less than in the ATLAS12 calculation. As a re- 
sult the line profile calculated with COSI appears slightly 
broader than the ATLAS12 calculation. However, the pro- 
files are practically identical when related to the same con- 
tinuum level. 

The same calculation over the wavelength range from 
200 to 1000 nm is shown in Fig.[3l In panel (a) the contin- 
uum irradiance spectrum calculated with ATLAS 12 (dot- 
ted line) and with COSI (solid line) are shown, where in 
the latter the hydrogen lines are included. Panel (b) shows 
the corresponding synthetic spectra including all spectral 
lines. To allow a better comparison the high resolution 
line spectra have been convolved with a 1-nm boxcar fil- 
ter. Generally, the calculations show a good agreement, 
and the overall shape of the continuum spectrum is well 
reproduced. The strongest discrepancies are between ^^350 
and 430 nm. The differences between the spectra are due 
to missing molecular lines in COSI and due to different 
photoionization cross sections as evident from panel (a). In 
Fig.Elpanel (c), the ratio between the COSI and ATLAS12 
line spectra is given for each 1-nm bin. Below 420 nm the 
ratio indicates some deviation while for longer wavelengths 
the ratio is very close to unity. 

To further illustrate the COSI to ATLAS12 compar- 
ison. Fig. [5] shows the same continuum calculation as in 
Fig.[3K, but over a narrower wavelength range in the UV. 
The discrepancy between 210 and 250 nm (panel a) can 
be explained by different photoionization cross sections 
which lead to a higher continuum flux in the COSI cal- 
culation than in the ATLAS12 calculation. We tested the 
effect of increasing the Mgi2 photoionization cross sec- 
tion by a factor of 2 in COSI. This leads to a much bet- 
ter agreement between both calculations but is inconsis- 
tent with Opacity Project data. The lower continuum cal- 
culated with COSI between 290 and 300 nm is a conse- 
quence of a resonance in the photoionization cross sec- 
tion of Mgi 3 level (cf. Fig.[T|). However, this resonance 
is not confirmed by the observations taken with the Solar 
Stellar Irradianc e Compari son Expe r iment (SOLSTICE) 



( Rottman et al.l . [1993 : Woods et al . 1993f l. as the mea- 



surements agree better with the ATLAS12 calculations. 
This indicates that there is still a need for improved 
atomic data. Nevertheless, from 260 and 280 nm (panel 
a) the continua give comparable results. In Fig. [6] the LTE 
line spectra calculated with COSI arc validated against 
the SOLSTICE measurements. The fairly good agree- 
ment indicates that the spectrum synthesis with COSI in 
LTE is calculated reasonably well. However, from ^310 
to 320 nm COSI leads to a considerable higher irradi- 
ance than ATLAS 12 and the SOLSTICE observation. The 
missing molecular lines of the CN-band explain the dis- 
crepancy between ~385 and 390 nm. 
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Fig. 4. Line opacities for tlie -wavelength range from 200 
to 210 nm for the depth point -where T5oo=1- The huge 
number of frequency points necessary to account for the 
detailed line opacities is reduced to 10 values by means 
of the ODFs. This allows the inclusion of the line opac- 
ities into the NLTE radiative transfer -with a reasonable 
computational effort. 



4. NLTE calculations 

4.1. NLTE-ODFs 
4.1.1. Theory 

For a complete calculation of the radiative transfer all 
relevant processes, e.g. the bound-bound, bound-free, and 
free- free processes have to be taken into account. Due to 
the huge number of bound-bound transitions the opacity 
of all spectral lines is considerable and must be taken into 
account in the NLTE radiative transfer calculation. In 
COSI only the line transitions bet-ween the explicit levels 
given m Table fO are calculated in NLTE. Thus, the hue 
opacity of a large number of lines is not directly included 
in the solution of the NLTE radiative transport. Indirectly 
"we account for them by applying opacity distribution func- 
tions (ODFs) of the line o pacities. The concept of ODFs 
has been f irst suggested b-y I Strom fc Kurucz ( 19661 ). for a 
review see lCarbon (11984 ." To our knowledge, this concept 
has not yet been used with an iterative scheme in a NLTE 
calculation. 



4.1.2. Implementation of NLTE-ODFs 

For the calculation of the ODFs the high resolution line 
opacities need to be calculated over the wavelength range 
under consideration, which in our case is from 90 to 
400 nm. The computing time would be prohibitive to ac- 
count for each single line in the NLTE solution and there- 
fore we describe all the line opacities by distribution func- 
tions over a certain wavelength range AA (in our case 
10 nm) as a function of line strength. This means that 
over AA and for each depth point the line opacities are 
sorted according to their strengths. Then, for a certain 



number of opacity bins (in our case 10) we calculate the 
mean opacity Kj for each bin j as given in Eq.[TJ 

mean(sort(Ki))j —>■ Rj (1) 

with i being the frequency index of the high resolution 
opacities. The bins j are then transfered to evenly spaced 
wavelength points within the wavelength interval under 
consideration yielding NLTE-ODFs with a 1 nm resolu- 
tion. Fig.m shows the high resolution line opacities (dot- 
ted line) from 200 to 210 nm at the depth where T5oo= 1. 
Also shown are the calculated ODF (dashed line) and the 
binned mean opacity (solid line) again for the same depth 
point. Thus, we include values of averaged line opacity, 
representing the high resolution line opacity, in the NLTE 
calculation. 

To avoid the excess of the chromospheric emission for 
lines that are calculated in LTE, the line emissivity for lay- 
ers outward the temperature minimum is set to the value 
calculated from the Planck function at the temperature 
minimum. We are aware that this approach clearly un- 
derestimates the emissivity for many upper chromosphere 
lines and it will be improved in the future through the in- 
clusion of more levels in the NLTE calculation and/or by 
an approximative treatment of the non-LTE effects, e.g. 
usi ng the m o dified excitation temperatures as introduced 
bv lSchmutd (Il99ll ). 

Although the NLTE-ODF distribution function is cal- 
culated from LTE level populations there is an indirect 
effect of the NLTE population through changes of the ion- 
ization equilibrium, i.e. all of the level populations that are 
treated in LTE with respect to the NLTE ground state 
population (cf. Table[XT]) are indirectly affected by the 
NLTE effects. Thus the inclusion of the ODFs changes 
the line strengths of the emergent spectrum, which in turn 
leads to a different set of ODFs. Therefore, self-consistent 
ODFs can only be calculated iteratively. We have investi- 
gated to what extent a set of successive iterations of ODFs 
influences the population values of the explicit NLTE lev- 
els. We also studied how many iterations are necessary 
to find a solution for which the ODFs are practically un- 
changed in a next iteration. The changes of the population 
numbers calculated without and with the first set of ODFs 
are up to a factor 30. The next iteration leads to changes of 
1.5, the second and third iterations show maximal changes 
of a factor of 1.002. As the population numbers are prac- 
tically unchanged between the second and third set of 
ODFs, we conclude that two iterations are sufficient to 
find a self-consistent distribution of ODFs (iHaberreiter . 
2006) . 

4.1.3. Results with NLTE-ODFs 

In this section we present NLTE calculati ons carried out 



with C OSI and using the F1999 model bv iFontenla et al 



(Il999l) along with the F2007 abundance values, and the 
NLTE-ODFs as described in the section above. In Fig. [7] 
the effect of the NLTE-ODFs on the emergent contin- 
uum spectrum is shown. The dotted line represents the 
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Fig. 5. As Fig. [3^ for the wavelength range 200 to 300 nm (panel a), and 300 to 500 nm (panel b). The differences are 
due to different photoionization cross sections implemented in COSf. The lines in panel (b) are the hydrogen lines 
longward of the Balmer jump. 
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Fig. 6. As Fig.[3lD plus a comparison with the SOLSTfCE measurements (dashed line). The synthetic spectra are 
convolved with f-nm boxcar filters to reproduce the resolution of the observation. The higher flux between 210 and 
235 nm can be explained by missing continuum and line opacity. The discrepancy between 385 and 390 nm is due to 
the CN-band not calculated with COSI. 



continuum spectrum calculated without NLTE-ODFs, the 
solid line shows the spectrum with NLTE-ODFs. The im- 
plementation of the NLTE-ODFs leads to a considerable 
decrease of the flux from ^100 to 200 nm. This is due 
to the fact that the NLTE-ODFs increase the opacity. 
Consequently, the height where r equals 1 is higher up in 
the solar atmosphere where the temperature is lower than 
in the case for the pure continuum opacity. Thus, the UV 
fl-ux is decreased and the elements, in particular the met- 
als, are less ionized. The increase of neutral atoms in turn 
leads to an increase of the continuum opacity and thus to 
a decrease of the continuum flux. It is important to note 
that the negative hydrogen opacity strongly depends on 
the electron density. As the ionization decreases, less neg- 
ative hydrogen is formed due to a decrease in the electron 
density in some layers, leading to a slight decrease of the 
continuum flux in the optical wavelength range (FiglTh). 
Note that COSI calculates the negative hydrogen in full 



NLTE, i.e. in a simultaneous solution together with all 
other NLTE levels as given in Table lA.l"] 

Fig. [5] shows the corresponding line spectrum calcu- 
lated in NLTE with COSI together with SOLSTICE mea- 
surements and ATLAS12 calculation. In order to com- 
pensate for the missing opacities, which is also partly 
present in the COSI LTE calculation, the Doppler broad- 
ening was increased from 1.5km/s to 15km/s in the wave- 
length range from 170 to 210nm, from 1.5km/s to 8km/s 
from 210 to 250 nm and 1.5km/s to 3 km/s from 250 to 
400 nm. For the wavelength range from 300 to 400 nm the 
COSI calculation is in reasonable agreement with the ob- 
servation. To reproduce the observation it is clear that 
the implementation of the NLTE-ODFs is important, be- 
cause without them the emergent flux would be consid- 
erably higher, as already shown in Fig. [71 Nevertheless, 
the remaining issues about the missing opacities need fur- 
ther investigation. For future calculations the inclusion of 
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Fig. 7. Comparison of the continuum calculations carried out with COSI in NLTE and using the F1999 model atmo- 
sphere structure to analyze the effect of the NLTE-ODFs on the calculated continuum. The implemented NLTE-ODFs 
for 100 to 400 nm lead to a decrease of the continuum flux, which in turn leads to less ionization. The decrease in 
electron density leads to a reduction of negative hydrogen opacity and therefore to a slight increase in the visible 
continuum (panel b). 
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Fig. 8. Comparison of the synthetic spectrum calculated with COSI in NLTE using F1999 as atmosphere structure 
with the SOLSTICE measurements (dashed line) and the ATLAS12 calculation. The synthetic spectra are convolved 
with a 1-nm boxcar to reproduce the resolution of the observation. 



molecular lines and additional continuum opacities such 5. Discussion 

as CH, OH and CN photodissociation are necessary. 



Finally, Fig. [5] shows the comparison of the line spec- 
trum calculated in NLTE with COSI together with 
SOLSTICE measurements for the wavelength range from 
100 to 400 nm. The overall agreement between the calcu- 
lation and the observation is good. However, the contin- 
uum around the Lyman- a line is considerably lower than 
the observation. One of the reasons for this is clearly the 
missing contribution of the quiet Sun network and also en- 
hanced network to the calculated spectrum. Also, the fact 
that our approach underestimates the emissivity of the 
chromospheric emission lines, can further explain the dis- 
crepancy below 150 nm. Nevertheless, the general agree- 
ment above 150 nm, in particular the continuum slope, 
clearly shows that the inclusion of the line opacities is 
essential for the calculation of the solar spectrum in the 
UV. 



In the present version of COSI we have calculated the 
NLTE-ODFs with line opacities in LTE with respect to the 
ground state of the ion. For all heights above the temper- 
ature minimum the emissivity is set to tL,\ = ^A.T^inKL^A 
with B being the Planck function and L the depth point 
index of the atmosphere structure. In F1999 the layers 
around Tmin are in near- LTE conditions which makes our 
approach suitable for this and similar stellar atmosphere 
models. However, if we would calculate an atmosphere 
with strong departures from LTE, as e.g. SRPM305, our 
assumptions would lead to a substantial overestimation 
of the opacity and at the same time an underestimation 
of the source function above the temperature minimum, 
which in turn would lead to a much lower UV continuum 
intensity than in the observat ions shown in the calcula- 
tions by iFontenla et al. ( 2007 ). In order to address such 
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Fig. 9. Comparison of the synthetic spectrum calculated 
with COSI in NLTE with the F1999 atmosphere struc- 
ture (solid line) together with the SOLSTICE observation 
at solar minimum on April 3, 1997 (dashed line). An in- 
creased line broadening has been applied to account for 
missing opacities. 



structures the NLTE effects have to be accounted for in 

mo re detail. 

Haberreiter et al. presented a comparison of the 



calculated Lyman-a line strength as a function of time 
over tw o solar cycles with i ts observed strength as com- 



piled by I Woods et al 



(200^. Their calculations are based 
on the same F1999 COSI model as discussed above. They 
found that the synthetic Lyman-a line is about a factor 
of 2 too weak. The discrepancy has been found in part to 
be due to ambipolar diffusion of protons and neutral hy- 
drogen as implemented in F1999 and the later atmosphere 
structure SRPM305 bv lFontenla et all (|2007t ) but not yet 
in the COSI code. Another important contribution which 
leads to an increased Lyman a line strength is the expan- 
sion of the area of active regions when a line is for med 
higher up in the solar atmosphere (|Sch511et al. ■ l2008h . 



6. Conclusions 

We presented the NLTE radiative transfer code COSI 
for the calculation of the spectral solar energy distri- 
bution and validated its performance in LTE against 
ATLAS 12 calculations and observed solar spectra. The 
results show that COSI calculates correctly the spectra in 
LTE. Furthermore, we have introduced the new concept 
of NLTE-ODFs. This new technique allows all line opaci- 
ties to be included into the NLTE radiative transfer. The 
fast convergence of the NLTE-ODFs does not significantly 
increase the computing time. We conclude that by inclu- 
ding NLTE-ODFs into the radiative transfer code COSI 
the calculations of the solar spectra are improved sub- 
stantially. It is important to note that this new concept is 
not limited to specific atoms, and can in principle also be 
used in an iterative scheme to derive stellar atmosphere 
structures. 
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Table A.l. Levels of the atomic model that are explicitly treated in NLTE by COSI. Given are the ion, the level and 
the atomic numbers AN; electron configuration and spectroscopic term designation, the ionization wavelength Aion 
calculated from the measured energies from the NIST database, where * denotes ionization to an excited state of the 
next higher ion; the ionization potential i^ion of the ground state; the energy E of the level, which is the weighted 
mean in the case of split levels; the statistical weight g of the model configuration. Fi nally t he las t column denotes 
the so urce from which the pho t oioniz ation cross se ctions were taken. H ere J88 refers to IJohnI (|l988h . MIH to lMihalad 
(UmI), seat to lSeaton etaP (|l994h and KOE to lKoester et al.l (|l985l ). Furthermore, Hyd stands for the hydrogenic 
treatment and OP means that the photoionization cross sections are from the OPACITY PROJECT and in case of 
iron from the IRON PROJECT. 
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12119.24 




33200.68 


2 


OP 


Nai 


7 


11 


[Ne]4d 


2d 


14486.66 




34548.74 


10 


OP 


Nai 


8 


11 


[Ne]4f 


2po 


14567.25 




34586.92 


14 


OP 


Nai 


9 


11 


[Ne]5p 


2po 


15601.91 




35042.03 


6 


OP 


Nai 


10 


11 


[Ne]6s 


^S 


19690 85 




36372 62 


2 


OP 


Nail 


1 


11 


[Ne] 


^s 






0.00 


1 




Me: I 


1 


12 


[Ne]3s2 


is 


1621.51 


61671.02 


0.00 


1 


OP 


Mgi 


2 


12 


[Ne]3s3p 


3po 


2512.26 




21877.30 


9 


OP 


Me: I 


3 


12 


[Ne]3s3p 


Ipo 


3755.50 




35051.26 


3 


OP 


Mg 11 


1 


12 


[Nc]3s 


2S 






0.00 


2 




All 


1 


13 


[Ne]3s23p 


2p° 


2077.19 


48278.20 


74.71 


6 


OP 


Alii 


1 


13 


[Ne]3s^ 
[Nc]3s23p2 


^S 






0.00 


1 




Si I 


1 


14 


3p 


1527.92 


65747.76 


149.68 


9 


OP 


Si I 


2 


14 


[Ne]3s23p2 


ID 


1686.36 




6298.85 


5 


OP 


Si I 


3 


14 


[Ne]3s23p2 


is 


1991.88 




15394.37 


1 


OP 


Si I 


4 


14 


[Nc]3s3p3 
[Nc]3s23p4s 
[Ne]3s^3p4s 
Nc 3s2('^S')3d 
[Ne]3s23p'' 




* 1328.47 




33326.04 


5 


OP 


Si I 


5 


14 


3po 


3875.02 




39799.50 


9 


OP 


Si I 


6 


14 


Ipo 


4062.78 




40991.88 


3 


OP 


Sill 


1 


14 


2po 






0.00 


6 




Si 


1 


16 


3p 


1199 85 


83559 10 


107.74 


9 


OP 


Si 


2 


16 


[Ne]3s23p'* 


Id 


1347.48 




9238.61 


5 


OP 


Si 


3 


16 


[Ne]3s23p'* 
[Ne]3s23p3 
[Ne]3s23p3 
[Ne]3s^3p^ 
[Ne]3s23p5 
[Ne]3s23p'' 


is 


1632.08 




22179.95 


1 


OP 


Sii 


1 


16 








0.00 


4 




Pi 


1 


15 


4go 


1182.30 


84580.83 


0.0 


4 


Hvd 

J.J.J' VJ. 


Pii 


1 


15 


3p 






0.0 


9 




CI I 


1 


17 


2po 


961.51 


104591.0 


0.0 


6 


Hvd 


CI II 


1 


17 


3p 






0.0 


9 




Ar I 


1 


18 


[Ne]3s23p'5 
[Ne]3s23p^ 


is 


786.72 


127109.80 


0.0 


1 


Hvd 

J.J.J' VJ. 


Am 


1 


18 


2po 






0.0 


6 




Ki 


1 


19 


[Ar]4s 


^S 


2855.56 


35009.81 


0.00 


2 


Hyd 


Ki 


2 


19 


Ar 4n 


2po 


4547.06 




13023.64 


6 


Hyd 


Ki 


3 


19 


[Ar] 5s 


2S 


7149.41 




21026.55 


2 


Hyd 


Ki 


4 


19 


[Ar]3d 


2d 


7419.51 




21535.60 


10 


Hyd 


Ki 


5 


19 


[Ar] 5p 


2po 


9709.85 




24713.90 


6 


Hyd 


Ki 


6 


19 


[Ar]4d 


2d 


13132.87 




27397.50 


10 


Hvd 


Ki 


7 


19 


[Ar]6s 


2S 


13225.32 




27450.71 


2 


Hyd 


Ki 


8 


19 


[Ar]4f 


2po 


14526.59 




28127.85 


14 


Hyd 


Ki 


9 


19 


[Ar]6p 


2po 


16648.27 




29004.90 


6 


Hyd 


Ki 


10 


19 


[Ar]5d 


2d 


20721.96 




30185.40 


10 


Hvd 


Kii 


1 


19 


Arl 


is 






0.00 


1 




Cai 


1 


20 


[Ar]4s^ 


is 


2027 60 


49305 95 


0.00 


1 


OP 


Cai 


2 


20 


[Ar]4s4p 


3po 


2936.66 




15263.09 


9 


OP 


Cai 


3 


20 


[Ar]3d4s 


lp)0 


3641.13 




21849.63 


5 


OP 


Call 


1 


20 


[Ar]4s 


2S 






0.00 


2 




Sci 


1 


21 


[Ar]3d4s2 


2D 


1896.81 


52922.0 


0.0 


10 


Hvd 

J.J.J' VJ. 


Sc II 


1 


21 


[Ar]3d4s 

[Ar]3d24s2 

[Ar]3d24s 


3D 






0.0 


15 




Til 


1 


22 


3p 


1831.05 


55010.0 


0.0 


21 


Hvd 


Tin 


1 


22 


4f 






0.0 


28 




Vi 


1 


23 


[Ar]3d34s2 
[Ar]3d^ 


4f 


1861.26 


54360.0 


316.484 


30 


Hyd 


Vii 


1 


23 


^D 






0.0 


25 





Table A.l. Continued. 



Haberreiter et al.: NLTE calculations of the solar atmosphere 



Ion 


Level 


AN 


Configuration 


Term 


Aion (A) 


Eion (cm 1) 


E (cm-i) 


9 


cross 


Cri 


1 


24 


[Ar]3d'Ms 


'S 


1832.32 


54575.6 


0.0 


7 


Hvd 


Cr II 


1 


24 


[Ar]3d^ 








0.0 


6 




Mm 


1 


25 


[Ar]3d'Ms2 




1667.80 


59959.40 


0.0 


6 


Hvd 


Mn II 


1 


25 


[Ar]3d'''4s 








0.0 


7 




Fei 


1 


26 


[Ar]3d'54s2 




1589.04 


63737.00 


402.96 


25 


OP 


Fei 


2 


26 


[Ar]3d7(4F)4s 


5f 


1789.67 




7457.82 


35 


OP 


Fei 


3 


26 


[Ar]3d^(4F)4s 


3f 


* 1873.48 




12408.10 


21 


OP 


Fei 


4 


26 


[Ar]3d7(4p)4s 
[Ar]3d'54s2 


5p 


2193.69 




17761.40 


15 


OP 


Fei 


5 


26 


3p 


2266.85 




19232.20 


9 


OP 


Fell 


1 


26 


[Ar]3d'5(5D)4s 
[Ar]3d'^4s^ 


6D 






0.00 


30 




Co I 


1 


27 


4f 


1616.96 


63430.0 


792.80 


28 


Hvd 


Con 


1 


27 


[Ar]3d8 


3F 






0.0 


21 




Nil 


1 


28 


[Ar]3dS4s2 


3f 


1676.27 


61600.0 


971.80 


21 


Hvd 


Ni II 


1 


28 


[Ar]3d9 
[Ar]3d^04s 


2d 






0.0 


10 




Cui 


1 


29 


2S 


1604.69 


62317.44 


0.0 


2 


Hvd 


Cuii 


1 


29 


[Ar]3di0 


is 






0.0 


1 




Zni 


1 


30 


[Ar]3di04s2 


is 


1319.80 


75769.33 


0.0 


1 


Hyd 


Znll 


1 


30 


[Ar]3di04s 


2S 






0.0 


2 
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